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Abstract 

Oscillations of atomic nuclei in crystals are considered in this paper. It is shown that elastic nuclei 
oscillations relatively electron envelops (inherent, I-oscillations) and waves of such oscillations can 
exist in crystals at adiabatic condition. The types and energy quantums of I-oscillations for different 
atoms are determined. In this connection the adiabatic crystal model is offered. Each atom in the 
adiabatic model is submitted as I-oscillator whose stationary oscillatory terms are shown as deep 
energy levels in crystals. The I-oscillations can be created at the expense of recombination energy 
of electrons and holes on electron- vibrational centers in semiconductors. They interact among 
themselves, with phonons and electrons, they influence on physical properties of crystals and 
crystal structures. The I-oscillations representing oscillations of atomic nuclei relatively electron 
system in crystals or molecules are the important physical reality. 

PACS numbers: 63.00.00, 72.2Q.-i, 78.20.-e 



1 



I. INTRODUCTION 



Electronic properties of crystals (and molecules) usually investigate by using of ap- 
proached solutions of Schrodinger equation. The important results of such investigations 
are properties of electronic system but the available information about energy quantum 
and types of nuclei oscillations are limited because of features of used adiabatic approxima- 
tions. An information about crystal oscillations usually receive from the dynamic analysis 
of traditional crystal models in which atoms are replaced by physical points with masses of 
atoms. Such traditional models contradict adiabatic theories because nucleus and electron 
envelop of every atom in this model is one particle and adiabatic nucleus oscillations (rel- 
atively electron environment) are impossible. Information about nuclei oscillations which 
was received on basis of traditional crystal models is incomplete. At the same time the 
role of nuclei oscillations in various physical effects, such as phonon drag of electrons and 
hyperconductivity - superconductivity at room and at higher temperatures, is important. 
Scientific understanding of such physical effects is inconvenient without enough information 
about nuclei oscillations and interaction of these oscillations with electrons and phonons. 
In this connection it is useful to consider the nuclei oscillations in crystals in accordance 
with adiabatic theory of solids, to construct the adiabatic crystal model in agreement with 
adiabatic theories, to estimate possible influence of nuclei oscillations on physical properties 
of crystals or molecules. Given paper is devoted to these questions. 



II. THEORETICAL PRECONDITIONS 



The modern theory of solids contain the representation about atomic structure of sub- 
stances and Schrodinger equation: 

m- = ™. (i, 

Hamiltonian (H) for a crystal include the kinetic energy operators for nuclei (T z ) and for 
electrons (T e ), and also crystal potential (V): 

H — T z + T e + V. (2) 

Crystal potential V usually include the operators of Coulomb interaction between elec- 
trons (V e ), between nuclei (V z ), and between electrons and nuclei (V ez ) :V — V e + V z + V ez . 
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Solution of Eq. (|l|) is the wave function of a crystal ty(r,R,t) where t - time, r and R 
designate degrees of freedom for electronic system and for nuclei system, accordingly. It is 
impossible to solve precisely the given equation dependent on huge amount of (more than 
10 23 ) variables because of known fundamental and technical difficulties. Theory searches 
for the reasonable approached solutions of the Eq. (|I|) which would allow to calculate and 
to predict noticed on experience physical magnitudes, and crystals properties. In this con- 
nection apply an approached methods. So, division of variables in the Eq. ([!]) would allow 
to simplify the problem and to consider motion of electron system irrespective of nuclei 
system. In such case these systems would submit to micro-canonical distribution, and the 
exchange of energy between systems would be impossible, as though they were divided by 
adiabatic membrane. For the first time P. Dirac has applied separation of electronic and 
nuclear variables |l[ and has presented wave function as ^(r, R, t) = <fi(r,t) ■ x(R,t)- He 
show that, supposing the energy conservation, it is possible to replace the Eq. (JJ) by two 
equations for wave functions 0(r, t) and x{Ri t), which describe movements of electronic and 
nuclear systems of a crystal: 

= ~ E |^ V ^ + {J dR - X*(Ri, t) ■ V(r h Rj) ■ X (Ri, t)}4>, (3) 

i 

i 1 j 

where m - electron mass, Mi - mass of I-th nuclei, V(ri, Ri) - crystal potential and H e (ri, Rj) 
- Hamiltonian of electron system. Coupled equations introduce the basis of the time- 

dependent self-consistent field (TDSCF) method. One can see from the Eqs. ©-© that the 
potential field, in which each system of particles move, depends on time and from results of 
averaging on coordinates of other system. In it the new quality consists which electrons get 
at the presence of nuclei. That is visible even on example of hydrogen atom, where proton 
and electron separately have coulomb fields, but in envelop of hydrogen atom the electron 
move in coulomb field of nucleus, and the nucleus move in a field of electronic shell, which 
have not coulomb field. The electronic envelops of atoms can unite among themselves, form- 
ing molecules and crystals but electrons outside of shells can not do that. So the electronic 
envelops represent new collective quality of electrons in atom but exact division of variables 
in the Eq. ([I]) cannot be carried out, and the problem about movement of nuclei system and 



3 



electron system in a crystal generally is self-consistent and nonadiabatic. Nevertheless adi- 
abatic approximations are used. Born-Oppenheimer adiabatic approximation [[J is related 
to the solution of the stationary Schrodinger equation 

= WV, (5) 

where W - energy of a crystal. Wave function of a crystal is product of electronic and nuclear 
wave functions: 

V(r,R) = (p(r,R)-$(R). (6) 
One can receive from the equation Eq. (|5]) the following two equations 

(T e + V) ¥ >(r,R) = E<p(r,R), (7) 



(T Z + E + A)$(R)=W$(R), (8) 
where E - energy of electron system and 

A = - J2( h2 / 2M i) [ <P*(r, R)V 2 Mr, R) ■ dr r (9) 
i 

- adiabatic potential. Equations (0, |8|) form the connected system because potential A in- 
cludes dependence of electronic wave function from change of nuclei coordinates and crystal 
potential V depends from r and R. In it the reasons of unadiabatic interaction of electrons 
system with nuclei system consist. Potential A describe stationary and, hence, oscillating 
process of energy exchange between systems of electrons and nuclei. However adiabatic 
conditions are provided with the certain accuracy, if potential A is small or is equal to zero. 
The Eqs. under condition of A = are known as Born-Oppenheimer adiabatic ap- 

proximation. The accuracy of the given approximation is determined by the contribution of 
potential A in energy of a crystal. The appropriate contribution of displacements, rotations 
and oscillations of nuclei system as a power series of small parameter rj = (m/M) 1 / 4 was 
calculated in Ref . || . In accordance with the energy of nuclei oscillations is proportional 
rf. In this connection sometimes assert, that the small magnitude of the relation (m/M) 
ostensibly is a condition of a validity for using of adiabatic approximation. This opinion 
is inexact, as the relation (m/M) is always small, but the validity of adiabatic approach 
depends on physical conditions and not always is justified. The wave function in Eqs. (^) 
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generally speaking differs from the exact solution and after her substitution in the Eqs. ([5]) 
the additional components turn out due to which the movements of nuclei are capable to 
cause electronic transitions and by that to break adiabatic condition. C. Herring has shown 
P|, that adiabatic approximation is correct, if 

/d 
d 3 r V *—Vj, (10) 

where E^ - energy of enabled electronic transition between conditions i and j, v - number 
of oscillation modes for nuclei system, dR v - characteristic displacement of nuclei system at 
frequency with number v. In case of crystals the systematic displacement and rotations of 
all system of nuclei can be excluded. It may to take into account only oscillations of nuclei 
at considering fixed crystal for estimation of adiabacity. Then in agreement with || the 
sufficient condition for validity of adiabatic approach on the only frequency with number v 
is smallness of nuclei oscillations energy in comparison with E^: 

Eij > hu u . (11) 

Believing that the conditions Eqs. (]T0|)-(]TT|) are satisfied, sometimes, one apply additional 
approximation in which wave function of a crystal (for estimation of potential A and for 
calculation of energy for electronic system E) is product of independent from each other 
wave functions of electronic and nuclear systems: 

*{r,R) = <p{r,Ro)-*{R), (12) 

where Rq - set of nuclei equilibrium positions. 

Thus, the stationary oscillations of nuclei system are possible at adiabatic approximation 
conditions. Adiabatic oscillations of nuclei system do not depend on oscillations of electronic 
system. They do not influence electronic system and do not change electronic structure of 
a crystal as a whole. From the point of view of adiabatic theories the crystal lattice is 
formed by system of electrons, and the nuclei are placed in potential minima of electronic 
system, where they can carry out oscillations concerning electronic system. The opportunity 
of nuclei oscillations becomes a reality when the necessary energy is transmitted to them. 
These conditions can be executed, for example, with the help of the electron-vibrational 
centers (EVC). Adiabatic oscillations of nuclei are important but for the present time they 
are a poorly investigated type of oscillations in crystals and molecules. 
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III. OSCILLATIONS OF ATOMIC NUCLEI 



Solutions of the Eqs. (0) and (H) m adiabatic approximation do not depend from each 
other. The analysis of nuclei system adiabatic oscillations may be carried out by analyzing 
the solutions of Eq. (|8|). This many particle problem about movements of nuclei system 
is possible to reduce, for example by Hartree method |J, to one particle problem about 
movement of one j-th nucleus in effective potential field V(Rj) depending from coordinates 
only of j-th nucleus: 

T j + V(R j )F j = W j F j , (13) 

where Tj - operator of kinetic energy and Fj - wave function, Wj - energy of stationary 
oscillations of j-th nucleus in a potential field V(Rj). This field is created by all electrons 
and nuclei of a crystal besides the j-th nucleus. The electronic environment of j-th atom 
brought basic contribution in V(Rj) . Contribution of electrons and nuclei from other atoms 
of a crystal is insignificant because of their symmetric and removed position concerning j-th 
atom nucleus. Adiabatic oscillations of nucleus occur on frequencies with number v (which, 
generally speaking, differ from characteristic frequencies of electronic system oscillations) 
near to a minimum V(Rj) conterminous with the center of atom electronic shell. They do not 
cause an electronic transitions and the electronic environment of each atom remains constant 
and motionless on frequencies with number v. It is possible to name these oscillations as 
inherent or I-oscillations, because of their properties are determined by inherent parameters 
of atom: mass and charge of a nucleus, potential V(Rj) near to the center of an electronic 
envelop. 

Accordingly each atom in adiabatic models of a crystals should be presented as inherent 
oscillator (I-oscillator) consisting from a nucleus and electronic shell, which is connected 
with each other by elastic force. In such case unadiabatic processes are transitions between 
stationary conditions of adiabatic model. The electron-vibrational transitions caused by 
presence of potential A in the Eq. (H) are a typical example of unadiabatic processes in 
which I-oscillations of the electron-vibrational centers, electrons and phonons strongly inter- 
act among themselves. Thus the I-oscillator energy levels (I-terms) of the centers are shown 
as deep energy levels. Nuclei oscillations relatively electrons shells of atoms in crystal occur 
in small area (~ 10 2 A) and in general they are necessary to research by quantum methods. 
Nevertheless in classical molecular dynamics the nuclei movements describe by Newtonian 
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equations of motion. If to study harmonic oscillations then it is suitable to use the known 
correspondence between results of quantum and classical theories. This correspondence 
consist in coinciding the transition frequency between adjacent quantum oscillatory levels 
of harmonic oscillator with his classical oscillation frequency. So energy spectrum of har- 
monic I-oscillations is possible to research by classical method. We used this possibility for 
description of oscillations in adiabatic crystal model when nuclei and electron envelops can 
move along straight axis U. Such diatomic model is shown in the top part of Fig. [I] where 
circles are electronic envelops and closed points in circles centers are nuclei of atoms. System 
equations of motion for this model may be written: 
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M ^ U 'ln = -h(U' ln -Ul), (14) 

M 2 ^U' 2n = -b 2 (U' 2n -Ul), (15) 
d 2 

m iQ^ u in = ~h{uL ~ u \n) ~ 9i(U' ln - U' 2n ) - g 2 {U' ln - Ul n _ x ), (16) 
d 2 

m2 Q^ U 2n = - h ^ U 2n - u' 2n ) - gx(U 2n - U ln ) - g 2 {U 2n - u' ln+1 ), (17) 

where Mi and M 2 - nuclei masses, mi and m 2 - an electronic envelops fictitious masses, U 
and U" - displacements of nuclei and electronic shells from equilibrium position , bi, b 2 , g\, 
g 2 - coefficients of elastic forces, t - time, n - elementary cell number. If to be interested 
in waves of harmonical oscillations in a linear atomic chain consisting of atoms which differ 
only by factors g\ and g 2 , it is necessary to put m\ = m\ = m , Mi = M 2 = M, b\ = b 2 = b. 
Let's accept for definiteness g\ > g 2 . Then the following dispersion relations are true: 

u J t 2 ^M = (D/2)±(D 2 /A-Fy/ 2 , (18) 

where D = [M(b + g x + g 2 - C) - m6](mM)~ 1 ; F = [b(g 1 + g 2 - B)](Mm)- 1 , B = -[g x + 
g 2 + 2gig 2 — C] 1 ^ 2 , C = 1 — 2sin(ag/2), q - wave vector, a - lattice spacing. The dispersive 
curves qualitatively are represented in Fig. [l] and contain known acoustic (A) and optical 
(O) branches, and also new branches (I) which concern to inherent oscillations representing 
oscillations of nuclei relatively electronic envelops. In this connection it is expedient to 
determine magnitudes of elementary quantums for I-oscillations. Nucleus of j-th atom (or 




FIG. 1: Dispersive branches of acoustical (A), optical (0) and inherent (I) oscillations for linear 
two- nuclear chain according to Eq. (|l8|), A fragment of such chain (cells with numbers n and n+1) 
is shown in the top part of the figure. 

ion) of crystal oscillate in effective field V(Rj). It is obvious that the basic contribution in 
V(Rj) introduce the electrons from envelop of considered j-th atom but contribution of nuclei 
and electronic envelops both environments of the next and more distant atoms is insignificant 
because of their symmetric and distant location relatively of j-th atom. This feature together 
with deduction of adiabatic theory allow to calculate frequencies of nucleus I-oscillations in 
a field of his motionless electronic shell. The potential function V(Rj) depends mainly 
on s-electrons density near the center of electronic envelop and in first approximation is 
spherical symmetric. Therefore small oscillations of nucleus near to the center of electronic 
envelop (a-type of I-oscillations) practically do not depend on polarization and have small 
unharmonicity. Quantums of I-oscillations in a crystal can be calculated when the function 
V(Rj) is known. So, elementary quantum of harmonic inherent oscillations in statistical 
model of Thomas- Fermi atom does not depend on nuclear number and is equal to 0.43 eV. 
This energy can be considered as majorization of elementary quantum for I-oscillations 
in multi-electronic atom. Function V(Rj) can be found for any atom from solving the 
Poisson equation with electrical charge density created by electrons of envelop near to her 
center. In atom of hydrogen the electronic density is defined by wave function: e|\l/i s | 2 , 
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where e - electron charge. The appropriate magnitude of quantum for harmonic inherent 
oscillations Uu\ = 0.519 eV. Un-harmonic correction to huj\ in the first and second orders 
of the perturbation theory in condition with oscillatory number v = is equal to +7 meV 
(2.7 percent), in condition with v = 1 correction is equal to +35 meV (4.5 percent) and 
in a condition with v = 2 correction is equal to +70 meV (5.4 percent). The similar 
calculations for helium give magnitude of quantum hu>2 = 0.402 eV. Electronic density near 
to the center of electron envelop in multi-electronic atoms is created basically by s - electrons. 
Quantums of inherent oscillations (a - type of I-oscillations) for atoms with numbers Z > 80 
were calculated taking into account the s-conditions and the shielding a nuclei by electrons. 
Quantums of a - type I-oscillations for atoms with Z > 2 may be determined with help of 
the following formula 

where H = 1.2 takes into account the contribution in electronic density from 2s- and 3s- 
conditions, £ = aZ 1 ^ 3 , magnitude a changes from 1 up to 1,15 at increase Z from 2 up 
to 80. Various types of I-oscillations in a crystal connected to various parts of atoms are 
possible. At increasing Z the electronic density at the center of atom envelop is considerably 
increased when emerge s-electrons of L -, M -, N - orbitals. Therefore it is expedient to 
define frequencies of I-oscillations for atom parts distinguished by the contents s-electrons. 
So, the joint oscillations of a nucleus and L - electrons (nucleus and two Is electrons) 
concerning other part of environment represent /3-type I-oscillations. Elementary quantum 
of /3-type I-oscillations may be calculated for atoms with Z > 2 with help of the Eq. (|T9"D 



at s — o{Z — 2) 1 / 3 and H = 0.2. It is possible to calculate elementary quantum of 7 - type 
I-oscillations representing joint oscillation of a nucleus, K - and L - orbitals concerning other 
part of an environment for atoms with Z > 8 with help of the Eq. fll9|) at s = o{Z — 8) 1 / 3 
and H = 0, 056. Thus it is possible to describe terms of I-oscillations a—, j3— and 7— types 
by the formula for harmonic oscillator 

?iw v = hu z (u+ 1/2), (20) 



where oscillatory quantum number v — 0, 1, 2, ... and hu z is defined by the Eq. (|19D . These 
oscillatory energy levels together with the origin of energy scale E = are energy levels of 
adiabatic crystal model and correspond to adiabatic (a— ,j3— ,7— type) nuclei oscillations 




FIG. 2: Calculated dependencies of the elementary quantum energies for I-oscillations of a- type 
(a), of /3-type (/?), and of 7-type (7) from atomic number Z are shown by open points. Experimental 
quantum energies of I-oscillations for several atoms are shown by closed points. 

concerning electron system in a crystal. The transitions between the energy I-oscillatory 
levels or on the energy I-oscillatory levels occurring with electron and phonon participation 
are non-adiabatic processes which may be studied with the help of offered adiabatic crystal 
model. The calculated quantum meanings depending on Z are shown in Fig. 2 by open 
points for a-type, /3-type, 7-type of I-oscillations. Experimental quantum meanings of I- 
oscillations for atoms: carbon, oxygen, sulfur, aluminium, iron, gadolinium, which were 
measured by various methods. They are shown in Fig. 2 by closed points. 
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IV. CREATION AND DESTRUCTION OF INHERENT OSCILLATIONS 



Inherent oscillations can be created at the expense of suitable energy source for example 
thermal, optical, energy recombination in semiconductors. The thermal way is represented 
improbable, as the Debye temperatures for inherent a-type oscillations (more than 2500 K) 
exceed temperatures of melting for many materials. The way for excitation of I-oscillations 
at the expense of energy allocated at recombination in semiconductors probably is quite 
effective and practically important when recombination occurs on the electron-vibrational 
centers (EVC). The EVC in semiconductors usually represent an association of impurity 
atom with vacancy and are characterized by strong interaction of electrons with phonons. 
In the recombination act on such centers participate average S > 1 phonons, where S 
- constant for electron-phonon interaction (Pekar-Huang-Ryse constant) which can reach 
several tens. The EVC equilibrium positions may be changed at recombination. That 
promotes to excitation of I-oscillations. Thus, the oscillations terms of EVC are shown 
as deep energy levels. The I-oscillation of atomic nucleus can formally be characterized 
by temperature T v = hu v /k , where k - Boltzmann constant. T v significantly exceed a 
temperature of a crystal lattice which are formed by electronic system. Or else, T v exceed 
temperature of electronic system. EVC disseminates this oscillatory energy in external 
environment with the greatest probability as crystal phonons, I-oscillations and waves of I- 
oscillations of the basic substance of a crystal. In result the recombination on EVC appears 
as the reason of I-oscillations and waves connected with phonons and capable influence on 
physical properties, of a crystal. Therefore it was possible to expect influence of EVC I- 
oscillations, and also I-oscillations of the basic substance atoms, on physical properties of 
semiconductor crystals and crystal structures. 

V. DESCRIPTION OF EXPERIMENTS 

In experiments were used flat polished semiconductor samples by thicknes 200 mkm con- 
taining EVC. It was GaP samples with impurities of aluminum or sulphur (~ 5 • 10 15 cm~ 3 ): 
GaP(Al) and GaP(S). Such impurities were chousen becouse atoms Al and S have masses 
appreciably exeed mass of atom Ga. This advantages emergence EVC and generation of I- 
oscillations. In experiments also was used quartz samples which where cut out from a single 
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crystal lengthways normal to axis c, and also samples of pyrolytic graphite and carbon nan- 
otube films on substrates which were fabricated by dispersion of graphite by electronic beam 
in vacuum ||. Carbon nanotube films represent regular structures from carbon nanotube by 
diameter about ten nm and in length about 0.1 mkm guided lengthways normal to a surface 
of a substrate. Silicon samples with impurities of phosphorus (~ 5 • 10 l5 cm~ 3 ) and oxygen 
(~ 10 18 cm -3 ): Si(P, O) were used also. In Si samples EVC was formed by oxygen atoms 
(A-centers: associations of oxygen atoms with vacancy, S = 5 for A-centers ). For perfor- 
mance of electrical measurements Au and Al contacts to samples were used with identical 
success. The contacts rendered by thermal dispersion of metals in vacuum. Temperature 
dependencies of electrical resistivity p{T) of samples containing EVC were measured within 
the range from 77 K to 700 K by a method Van-der-Pau and by two sonde method. The 
experimental temperature dependencies of resistivity in various sites are described by the 
formula p{T) = exp(E a /2kT). With the help of this formula were determined the activation 
energies (E a ) on various sites of experimental dependencies p(T). 

Infra-red reflection (IR) spectra connected with EVC were measured in the optical range 
from 15 mkm (83 meV) to 2 mkm (620 meV) at 300K. The angle between IR beam and flat 
surface of sample was equal to 45°. IR reflection spectra was analyzed on the basis of theory 
for reflection of classical charged harmonic oscillator J7|, IR reflectivity of the oscillator 

E - ^~ 1 ) 2 + fc2 (21) 

is defined by a parameter of refraction (n) and absorption coefficient (k). They depend on 
oscillator frequency (Q), frequency of crystal lattice oscillations (ujl), coefficient of damping 
(9) , optical permittivity (e op t) and optical frequency (uj) : 

" 2 - fc2 = £ ^ + " V-%-^ 2 ' 2nk = ul W^WT^- (22) 

The greatest meanings of are located between Q and u>l. Contour of the spectrum 
depends from e opt and 9. According to the quantum theory || [H| the experimental spectra 
of reflection was decomposed on components every of which fits spectrum for one oscillator 
by selecting parameters Q, u>l, e opt and 9. Thus the specified parameters were defined. 
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1000/T (K 1 ) 

FIG. 3: Temperature dependencies of resistivity (p) for undoped single crystal GaP (A), for GaP 
doped by aluminium (B) and for GaP doped by sulphur (C) are shown by open circles. Straight 
lines show tangents to the linear curves sites. 

VI. EXPERIMENTAL RESULTS AND DISCUSSION 



A. Temperature dependence of resistivity 

Typical temperature dependencies of electrical resistivity p(T) for GaP samples are shown 
in Fig. [| Magnitudes po were selected for each curve that the curves were conveniently 
arranged in the Fig. |3|. The sites of curves with constant inclinations are pointed out by 
pieces of line tangent to the curves. The inclination of lines tangents to the curves relatively 
axes of coordinates defines magnitudes E a . Magnitudes E a , measured at T < 330 K, are 
brought in Table I and basically correspond to known energies of phonons in GaP |JTT . 



Energy of phonons are introduced in the central column in Table I. Other energies in 
Table I supposedly may be connected with I-oscillations {(3- or 7-types) of Al or S and 
possible with combinations of such oscillations with phonons. 
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TABLE I: Activation energies for GaP(Al) and GaP(S) at T < 330 K. 



Activation energies E a (meV) for 
GaP(Al) samples with numbers: 


Phonons a 
in GaP 

(meV) 


Activation energies E a (meV) for 
GaP(S) samples with numbers: 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


8.1 






7.2 


8.0 


- 












15.0 






14.8 


14.5 


14.25 


15.0 


14.8 






14.3 


24.5 








25.0 


24.42 


24.0 




24.6 








35.0 
















28.0 


28.0 


45.0 








44.6 


44.75 




42.0 












49.0 


48.0 




47.00 


47.0 










75.0 




83.0 












70.0 







a Reference [11] 



The magnitudes E a for GaP (Al) and GaP (S) measured at T > 330 K are included in 
Table II. The Table II also contain energies calculated with help of Eq. ( |20| ) for inherent 
oscillations of Aluminum and Sulfur atoms with different v. The energies was calculated 
with hui3 = 0.283 eV for Aluminum and %U\q = 0.301 eV for Sulfur which were determined 
in section III and pointed out in Fig. ^| 

One can see from Table II that activation energies of samples with each type of impurities 
can be divided into two groups which corresponds to two columns dextral. One group 
consists of activation energies which are described by the Eq. (^(]). These activation energies 
are equal to energy of impurity atoms I-oscillations and correspond to transitions from 
oscillatory energy levels with v — 0, 1, 2, ■ • • in the minimum of oscillatory potential 
where oscillatory energy E = 0. Such transitions for free quantum harmonic oscillator are 
forbidden but they are possible due to interaction of I-oscillator with crystal. Consequently, 
inherent oscillators of impurity atoms show duality of properties that can be explained by 
their interaction with the environment. Other group of energies in Table II consists of ones 
multiple hu>i3 for Al and Hujiq for S. These groups of energies are also connected with inherent 
oscillations of impurity atoms and correspond to transitions between oscillator energy levels 
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TABLE II: Activation energies (eV) for GaP(Al) and GaP(S) at T > 330 K. 



E a for GaP(Al) samples 
with numbers: 


Calculated on 
Eqs. © 


Multiple 


1 


2 


3 


4 


5 






0.14 


0.14 


0.14 


0.138 


0.137 


0.142 (z/ = 0) 




0.28 


0.29 


0.29 


0.28 


0.28 


- 




0.42 


0.42 


0.43 


0.42 


0.43 


0.425 (u = 1) 






0.57 




0.56 


0.58 






0.71 




0.72 






0.707 = 2) 










0.85 










0.97 








0.991 (v = 3) 




1.10 








1.11 






E a for GaP(S) samples 
with numbers: 




Multiple 


6 


7 


8 


9 


10 






0.15 


0.15 


0.15 


0.15 


0.15 


0.151 (y = 0) 




0.30 


0.29 


0.30 


0.30 


0.31 






0.60 




0.60 


0.61 
















0.74 


0.753 = 2) 








0.92 










1.03 










1.050 (y = 3) 





with different v. One can see that huis = 0.283 eV is the same in both groups of energies 
for Al impurity, and Huiq = 0.301 eV is the same in both groups of energies for S impurity. 
Consequently both groups of energies refer to the same type of centers that show classical 
and quantum properties (duality of properties). 

The temperature dependencies p{T) for pyrolytic graphite (PG) samples also are char- 
acterized by several activation energies. The magnitudes E a for samples PG are submitted 
in Table III. They were measured at directions of a current parallel to nuclear planes (L) 
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TABLE III: Activation energies (eV) of resistivity for graphite samples at current directions parallel 
(L) or perpendicular (T) to atomic planes. 



Activation energies at L and T current directions 
and supposed type of oscillations (Type) 


Calculated 
on Eqs. @ 


Multiple 


L 


T 


L 


T 


L 


T 


T 


L 


Type 


- 


- 


- 


0.004 


- 


- 


- 


- 


0.003 


0.004 


TA 


- 


- 


- 


- 


- 


- 


- 


0.006 


0.006 


- 


TO 


- 


- 


- 


0.029 


0.021 


- 


- 


- 


0.024 


- 


- 


- 


- 


0.040 


0.040 


- 


- 


- 


- 


- 


- 


TA 


- 


- 


- 


0.044 


- 


- 


0.043 


- 


- 


0.043 


TA 


- 


- 














0.055 




TO 


- 


- 


0.082 




0.082 


0.084 




0.090 






LA 


- 


- 


0.125 




0.125 


0.125 


0.125 


0.125 


0.125 


0.125 


la 


0.125 {v = 0) 




0.250 


0.250 


0.247 


0.245 




0.250 


0.252 


0.250 


la 




hujQ 












0.374 






la 


0.375 {v) = 1 




0.500 




0.507 




0.500 








la 


















0.621 


0.628 


la 


0.625 (y = 2) 










0.759 










la 







and perpendicular to nuclear planes (T) of graphite. Also in Table III the energies of I- 
oscillations for carbon are brought. They calculated with help of the Eq. (P0| ) for various 
v. These energies was calculated with Tiloq = 0.25 eV for carbon which was determined in 
section III and pointed out in Fig. ^| 

Top lines of Table III contain of energies which may be identified with phonons of graphite 
in the certain points of Brillouin zone appropriate to greater density of phonons. It is possible 
to explain energies 0.021 eV, 0.024 eV and 0.029 eV in Table III by a combination several 
acoustical phonons. 

Activation energies located in the bottom of Table III can be divided into two groups. The 



first group consist of energies which close to energies calculated with help of Eq. (p0|) . Second 
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TABLE IV: Activation energies (eV) for carbon nanotube films on quartz and fluorite substrates. 



Nanotube films on quartz 
substrates with numbers: 


Nanotube films on fluorite 
substrates with numbers: 


Calculated on 
Eq. (|20|) 


Multiple 


1 


2 


3 


4 


5 


6 


7 


8 


- 


- 


0.013 


0.012 


0.013 


- 


- 


- 


- 


- 


- 


- 


0.017 


0.017 


0.017 


0.017 


- 


- 


- 


0.016 


- 


- 


0.030 


0.032 


0.030 






0.029 




0.035 


- 


- 


0.040 


0.050 


0.040 


0.037 


0.043 


0.044 




0.056 


- 


- 


0.095 






0.067 


0.095 


0.086 










0.127 


0.125 


0.125 




0.128 


0.121 


0.125 


0.125 


0.125 (v = 0) 




0.027 


0.250 


0.250 


0.250 


0.248 


0.250 


0.250 


0.250 








0.375 










0.375 




0.375 (v = 1) 










0.500 


0.495 




0.500 


0.500 


















0.875 




0.875 {y = 3) 





group consist of energies multiple Eq. The activation energies may be explained as expense 
of energy in the acts of charge carriers generation at transition EVC from a condition with v 
= 0,1 and 2 in minimum of potential V(Rj) where oscillatory energy E = 0. Others energies 
are multiple E and may be explained by transitions EVC between oscillatory condition with 
various v. These results specify duality properties of carbon I-oscillator in graphite similar 
duality of properties for I-oscillators in GaP(Al) and GaP(S). 

Experimental activation energies for carbon nanotube films on substrates are given in 
Table IV together with the calculated with help of the Eq. (|20"D data with subject to I- 
oscillator quantum of carbon atom {uiq = 0.25 eV). 

In the top lines of Table IV are located the energies measured at T < 300 K which can be 
identified with phonons and their combinations. The bottom lines of Table IV contain the 
activation energies which were measured at T > 300 K. This energies can be divided into 
two groups. One group consist of the energies which close to calculated energies for carbon 
I-oscillations with help of Eq. (|20|). It gives the certain basis to connect them to transitions 
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of carbon I-oscillators from condition with v = 0, 1, 2, 3 in a condition with oscillatory energy 
E = 0. Other group consist of energies multiple Twoq = 0.25 eV which can be connected with 
transitions of the same I-oscillators between oscillatory condition differing by v on 1 or 2. 

The results of measurements and analysis of activation energies in samples Si(P, O) are 
qualitatively identical to the stated results, but the activation energies at T > 300 K close 
to energies of oxygen I-oscillations (in A-centers) at v = 0, 1, 2, 3, 4. 

It is impossible to explain the experimental dependencies p{T) with scattering of charge 
carriers by phonons at T < 300 K since the scattering is capable to create an opposite 
effect to the observed effect of reducing resistance when increasing the temperature. We 
connect experimental activation energies with creation of free charge carriers at the expense 
of oscillations energy of EVC for example EVC created by atoms Al or S in GaP, by atoms 
C in graphite or in carbon nanotube films, and by atoms O (A-centers) in silicon. The 
coincidence of calculated and experimental activation energies at T > 300 K allows us 
to connect them with I-oscillatory terms of impurity atoms forming EVC. These terms are 
shown as deep energy levels. 



B. Infrared spectra 

Typical spectrum of reflectivity change (dR) in GaP(Al) is shown in Fig. §by curve A. The 
dR caused by introduction of impurity atoms (Al) in GaP in concentration ~ 5 ■ 10 15 cm -3 . 
According to the theory ||, |l(J the given spectrum was decomposed into components, which 
in Fig. 4 are denoted as B, C, D, E. The sum of the calculated spectra B, C, D and E, 
coincides with the experimental spectrum A if energies Tilou (i = 1, 2, 3, 4 ) coincide with 
energies of Al a -type I-oscillators: 0.5huJi 3 , hu>i 3 , l.hTiuJi?,, 2^13; hu± 3 = 0.283 eV. It is 
possible to see from these data that Al I-oscillators show duality of properties because in 
optical transitions the energy level of "zero oscillations" (0.5^13) is shown. It is possible to 
connect energy HQ = 61 meV with 7- type of atom Al I-oscillations. The attenuation (6) in 
I-oscillators is great: (0/Q) = 0.09. The best consent between experimental and calculated 
spectra is reached if e opt = 2 though for GaP e opt = 8.457 Hl2| . It may be the e opt = 2 
is necessary to carry not to a crystal GaP, but only to EVC. Thus the section of photon 
capture by EVC can depend on wavelength of phonons cooperating with EVC. The similar 
results were received at research of EVC reflectivity spectra formed by impurity atoms of 
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FIG. 4: Spectrum of reflectivity change (dR) which caused by impurity Al in GaP (A) and compo- 
nents of the spectrum connected with different I-oscillations of Al (B, C, D, E). Oscillator energy 
(HO,) and energies of Al a-type I-oscillations in different oscillations states hun (i = 1 5 2,3, 4) are 
shown by arrows. 

sulfur for which hu\§ = 0.301 eV. 

IR reflectivity spectrum of carbon nanotube film on molybdenum substrate (curve A) is 
shown in Fig. |5|. Spectrum A was decomposed on reflection bands of two oscillators with 
energies: hQ = 60 meV, hu^i = 0.25eV and hu^ 2 = 0.375eV. The first of the energies 
is equal to elementary quantum of carbon /3-type I-oscillations, second energy is equal to 
quantum of carbon atom a-type I-oscillations (hue) and last of the energies is equal to 
carbon a-type I-oscillations at v = 1. Thus, the experimental data about interaction a- and 
/9-types of I-oscillations and their influence on IR reflection spectra were received. 

The typical reflectivity spectrum of carbon nanotube film on a copper substrate contains 
one band which is well described by the calculated reflectivity spectrum of oscillator with 
parameters hQ = 60 meV and hu^ = 0, 25 eV which are accordingly equal to elementary 
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FIG. 5: Reflectivity spectrum of carbon nanotube film on molybdenum substrate (A) and calculated 
(at huiL2 = 0.375 eV, v = 1) reflectivity spectrum for carbon a-type I-oscillator (B). Oscillator 
energy and energies of carbon a-type I-oscillations are shown by arrows. 

quantums (3- and a-type I-oscillations of carbon atom. 

Reflectivity spectrum of single-crystal quartz is shown in Fig. 6 (curve A). Curve B 
in Fig. P represents calculated reflectivity spectrum of oxygen a-type I-oscillator with the 
following parameters: HQ = 0.11 eV; (ujl/Q) = 0.25 and e op t = 1.2, {0/Vt) = 0.011. One can 
see that energy 0.11 eV is equal to hui$/2 at v = in Eq. (p0|). Frequency ujl can be identified 
with frequency of a phonon. The satisfactory consent of spectra A and B in a Fig. 6 confirms 
the electron-vibrational nature of the given reflection band and allows definitely to connect 
her with excitation of a-type I-oscillations of oxygen atom at v — 0. Some distinction of 
spectra A and B is explained by the contribution to reflection a components with ul = 2Q 
and ujl = 3fi. 

Typical experimental reflectivity spectra of carbon nanotube film on quartz substrate 
are submitted in Fig. [7] (curves A, B, C). Spectrum A is measured at once after cultivation 
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FIG. 6: Reflectivity spectrum of quartz (A) and calculated reflectivity spectrum of oxygen a-type 
I-oscillator (B). 

of a film. Spectra B and C were measured 2 and 6 months later accordingly. Between 
measurements the film was stored at room conditions. The IR reflection factor (R w ) increase 
considerably during a storage of samples. The similar changes of reflection spectrum for 
carbon nanotube film occur at her heating in vacuum. The changes of a spectrum occurring 
during two months of storage in room conditions can be achieved at heating of a samples in 
vacuum at T ~ 500 K approximately within one hour. It allows to connect changes of the 
spectrum with changes of film structure. The experimental spectra contain a minimum in 
that area where the maximum of quartz reflection band is located. This reflection maximum 
is characteristic for impurity oxygen in silicon where it is displaced from 0.13 eV to 0.11 eV 
at increase of oxygen concentration and formation of quartz disseminations. This reflection 
maximum is characteristic also for silicon oxides, and for quartz. The reflection maximum 
definitely can be connected to a-type I-oscillations of oxygen atom at v — 0. Occurrence 
of reflection minimum close to 0.11 eV caused by carbon nanotube film on quartz substrate 
cannot be explained by anti-reflection coating action of a film because of the reflection 
minimum is observed at various thickness of a film and substrate. It is possible to explain 
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FIG. 7: Experimental reflectivity spectra of carbon nanotube film on quartz substrate after film 
fabrication (A), two month later fabrication (B), six months later fabrication (C) and reflectivity 
spectra of oxygen a-type I-oscillator (D, E, F), calculated at meanings of e op t equal to 2, 20, 250, 
accordingly. 

this minimum by interaction between I-oscillators of oxygen atoms in quartz substrate and 
carbon atoms in nanotube film. Curves D, E, F in Fig. |7| represent spectra of oscillator 
reflectivity calculated at parameters HQ = 0.11 eV and Tiujl = 1.25 eV with e opt equal to 
2, 20 and 250 accordingly. The first from specified energies coincides with energy of a-type 
I-oscillations for oxygen atom at v = in Eq. ( ^0|) (hug/2 = 0.11 eV) but the second energy 
is multiple to elementary quantum of a-type I-oscillations for carbon atom (hue = 0,25 eV). 
The reflectivity spectra of nanotube films on quartz after thermal treatment are satisfactorily 
described by oscillator spectrum at e opt = 10 3 — 10 4 but characteristic energy HQ = 0.11 eV 
and huji = 1.25 eV remain constant. It allows to explain the given experimental spectra by 
interaction with each other of oxygen and carbon I-oscillations at participation of phonons. 

In silicon samples the IR reflectivity spectra, optical transmission and photoconductivity 
at electron-vibrational transitions on A - centers have phonon structure and are described 
by line Pekar- Huang- Rhyse spectra JT3|, pi}] . Typical photoconductivity (curve A) and ex- 
perimental spectra of IR transmission (curve B) for containing A-centers (10 14 cm~ 3 ) Si 
single-crystal are introduced in Fig. [8|. Every of spectra lines correspond to participation 
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FIG. 8: Photoconductivity spectrum (A) and spectrum of IR transmission (B) measured at 80 K 
and connected with electron transitions on A-centers in silicon sample. Number of phonons which 
participate in electron- vibrational transitions is denoted as p. 

of p phonons. The energies of electron transitions without phonons participation (p = 0) 
determined on basis of the theory [|IB|, [TJ[] coincide with energies of a- type oxygen atom 
I-oscillations counted from top of valence energy zone in silicon. These spectra are measured 
in polarized IR radiation when the IR electrical vector was directed as normal to [100] in 
a crystal. The spectra contain extremes which differ from each other on energies multiple 
energy of optical phonon in Si (~ 55 meV). Height of spectral lines follows dependence 
(S p /p\), constant of electron-phonon interaction (Pekar-Huang-Rhyse constant) S ~ 5. En- 
ergy of electron transitions without phonons participation (p = 0) is close to energy £^+0.33 
eV (where E v - top of a valence energy zone) which correspond to v = 1 for oxygen a- type 
I-oscillator. At others polarizations of IR radiation it is possible to allocate spectra with 
participation others (optical and acoustic) phonons with the same meanings of constant 
S and energies of electron transition with p = conterminous with I-oscillatory terms of 
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oxygen atom. Curve B in Fig. || represents typical for EVC spectrum of negative photo- 
conductivity. Curve B in Fig. |^ represents spectrum of IR transmition. One can see that 
the increase of IR transmission (that is reduction of absorption) corresponds to photocon- 
ductivity increase. Thus terms of oxygen I-oscillations in A-center is similar to other atoms 
I-oscillations in structure EVC. The I-oscillation terms are shown as deep energy levels in 
semiconductors. 



C. Discussion 

Inherent oscillations of atomic nuclei relatively electron system (relatively electron en- 
velops of atoms) in crystals and molecules exist that one can see from results of experimen- 
tal and theoretical researches of different authors. The equations system Eqs. (|I4|)-(|T7|) at 
n = describes oscillations of one crystal cell that is oscillations of appropriate molecule. 
In correspondence with Eq. ( |T8D the oscillatory spectrum of a two-nuclear molecule contains 
frequencies of nuclei I-oscillations when electron envelops are motionless (on I-oscillations 
frequencies) and nuclei move relatively each other with inverse phase. Such oscillations re- 
ally are available in spectra of two-nuclear molecules. F. Vilesov and M. Akopiyan with 
help of optical mass spectrometry method investigated the ionization spectra for hundreds 



of molecules in various oscillatory conditions [15|]. In particular they defined the energies 
of transitions from the adjacent oscillatory conditions in molecule of oxygen (O2) whose 
differences are close to quantum of oxygen a-type I-oscillations (0.22 eV). That coordinates 
with I-oscillations presence in two-nuclear molecules. 



G. Pastore and E. Smargiassi have applied in Ref. |16] Car-Parrinello [17] and Born 



Oppenheimer [|18| molecular dynamics for calculation model system consisted of eight sil- 
icon atoms forming the periodic diamond lattice. Car-Parrinello dynamic basically takes 
into account conditions with various meanings of oscillatory quantum number u, but Born- 
Oppenheimer dynamic corresponds to the minimal oscillatory energy of a nucleus [y = 0). 
The comparison of application results for both dynamics allows to define transition ener- 
gies between conditions with various meanings v. In the Car-Parrinello dynamic "fictitious 
masses" about some hundreds of electron mass are used. They are inertia parameters as- 
signed to every orbital degrees of electron freedom at a nucleus motion. It is obvious that 
fictitious mass of electron envelop in every atom sufficiently exceed mass of nucleus. There- 
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fore at movement of a nucleus the electronic envelop remain motionless and changeless. It 
corresponds to inverted adiabatic approximation (inverted in relation to well-known adia- 
batic approximations) when nuclei move (hot nuclei) but the electrons "do not heat up" 
systematically in the presence of the hot nuclei (cold electrons). It corresponds to the 
representation about massive and motionless electronic envelop (on frequencies of nucleus 
I-oscillations) which we used for definition energy quantums of I-oscillations. 

In agreement with Ref. |16| the density of electronic conditions in considered atomic 
model system reaches the maximum near 1.3 eV that close to width of the forbidden zone in 
Si single crystal (E g = 1.16 eV at temperature K). The minimal energy in the calculated 
electron spectrum (0.6 eV) exceeds the maximal phonons energy (= 55- 10~ 3 eV) that satisfies 
inequalities Eqs. ([T0|)-([TTD and does improbable systematical flow of energy between nuclei 
and electrons. The applicability of adiabatic approximation is thus reasonable. In [It] also 
are calculated: kinetic energy of the electrons (T e ), electronic energy (V e ), kinetic energy of 
nuclei (T z ), conserved energy (E cons = T e +T z +V e ), physical total energy (E phys = E cons -T e ). 
Energies E cons and E p h ys are strictly constant. Energies V e and T e change on 0.43 eV 
and 0.81 • 10~ 3 eV accordingly in opposite phase to each other with the common period 
T\ = 2.15 ■ 10~ 14 s which correspond to energy E\ = (27rh/ri) = 0.194 eV. The specified 
energy 0.43 eV can be connected to transitions of silicon I-oscillator between condition in 
minimum of potential V(Rj) which coincide with oscillatory energy of nucleus E = and 
condition with v = 1 that is calculated with help of Eq. (^) for silicon atom with meaning 
huuu = 0.29 eV. If it is correct then calculations in |16| proves observable on experience 
the duality properties of I-oscillators. In accordance with [|l6l the electron system acts on a 
nucleus by force which in Car-Parrinello approach (Ecp) and in Born-Oppenheimer approach 
{Fbo) oscillate in phase with the common period t 2 = 4.32 ■ 10~ 13 s which correspond to 
energy E 2 = {2tiTi/t2) = 9.67 • 10~ 3 eV. Besides, the difference (Fcp — Fbo) oscillate with 
two periods which correspond to energies = 9.672 ■ 10 -3 eV and = 0.2998 eV. Last of these 
energies is near to energy of I-oscillator transitions for silicon atom between the adjacent 
oscillatory conditions with v = (Born-Oppenheimer approach) and v = 1 (Car-Parrinello 
approach) . 

One can see from these results that in adiabatic approximation the complete crystal en- 
ergy is constant but there is the oscillating process of energy exchange between system of 
electrons and system of nuclei. Therefore the adiabatic conditions Egs. (|TU|)-(|TTD concerning 
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systems of electrons and nuclei is correct only for average energy on an time interval exceed- 
ing the periods T\ and r^- The I-oscillations of atomic nuclei carry out important role in 
oscillating process of energy exchange between systems of electrons and nuclei in crystals. 
These I-oscillations give raise to such phenomena as phonon drag of electrons [T^], thermal 
superconductivity, and hypercoductivity at very high temperatures. 

VII. CONCLUSION 

Electronic envelop of any atom represents the new collective quality of electron system 
in comparison with quality of free electrons. The electronic envelops of different atoms can 
unite with each other and form crystal lattices due to that the nuclei of atoms get an oppor- 
tunity to carry out adiabatic inherent oscillations (I-oscillations) relatively electronic system 
of a crystal. In adiabatic approach there is no systematic energy flow from electronic system 
to nuclei system or in opposite direction. However there is the oscillating process of energy 
exchange between system of nuclei and system of electrons which is submitted by inherent 
oscillations in particular. These I-oscillations are a nucleus oscillations relatively motionless 
(on I-oscillation frequencies) nucleus environment in a crystal (a-type of I-oscillations). It 
can be also oscillations of a nucleus together with K-electrons (/3-type of I-oscillations) or 
nucleus together with K- and L-electrons (7-type of I-oscillations) relatively motionless (on 
frequency of I-oscillations) environment. Each atom in adiabatic model of crystal is sub- 
mitted as appropriate I-oscillator. Non-adiabatic processes can be considered as transitions 
between stationary terms of I-oscillators. 

Stationary I-oscillations and also waves of I-oscillations can exist in crystals. The energy 
spectrum of I-oscillations and waves of I-oscillations can be described by the formula for 
harmonic oscillator (at neglecting by non-harmonicity) . The elementary quantum of har- 
monic I-oscillation depend on nuclear number, and for a-type I-oscillations they considerably 
exceed energy of phonons. 

The effective creation of I-oscillations and waves of such oscillations can be carried out at 
the expense of recombination energy of electrons and holes on electron-vibrational centers 
in crystals. I-oscillations energy terms of the electron-vibrational centers are shown as deep 
energy levels in semiconductors. 

The I-oscillations influence on physical properties of crystals and crystal structures. Thus 
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the inherent oscillations representing oscillations of atomic nuclei relatively electronic system 
in crystals or molecules are the important physical reality. 
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